LARGE TIME DECAY AND GROWTH FOR SOLUTIONS OF A 
VISCOUS BOUSSINESQ SYSTEM 
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Abstract. In this paper we analyze the decay and the growth for large time of weak 
and strong solutions to the three-dimensional viscous Boussinesq system. We show that 
generic solutions blow up as t — > oo in the sense that the energy and the L p -norms of the 
velocity field grow to infinity for large time for 1 < p < 3. In the case of strong solutions 
we provide sharp estimates both from above and from below and explicit asymptotic 
profiles. We also show that solutions arising from (uo,9o) with zero- mean for the initial 
temperature 80 have a special behavior as or t tends to infinity: contrarily to the 
generic case, their energy dissipates to zero for large time. 



1. Introduction 

In this paper we address the problem of the heat transfer inside viscous incompressible 
flows in the whole space R 3 . Accordingly with the Boussinesq approximation, we neglect 
the variations of the density in the continuity equation and the local heat source due to 
the viscous dissipation. We rather take into account the variations of the temperature by 
putting an additional vertical buoyancy force term in the equation of the fluid motion. 

This leads us to the Cauchy problem for the Boussinesq system 

( d t 6 + u-V6 = kA8 
dtu + u • Vn + Vp = vAu + fiOes 
V-ii = 

u\t=o = Uq, 6\t=o = Go- 



(1.1) 



x S R 3 , t G 



Here u: R 3 x R+ -> R 3 is the velocity field. The scalar fields p: R 3 x R+ ->• R and 9 : M 3 x 
M + —J- M denote respectively the pressure and the temperature of the fluid. Moreover, 
es = (0,0, 1), and (3 S R is a physical constant. For the decay questions that we address 
in this paper, it will be important to have strictly positive viscosities in both equations: 
v, k > 0. By rescaling the unknowns, we can and do assume, without loss of generality, 
that v = 1 and (5 = 1. To simplify the notation, from now on we take the thermal diffusion 
coefficient k > such that k = 1. 
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Like for the Navier-Stokes equations, obtained as a particular case from (jl.ip putting 
8 = 0, weak solutions to (II. ip do exist, but their uniqueness is not known. The global 
existence of weak solutions, or strong solutions in the case of small data has been studied 
by several authors. See, e.g. pQ, [11], [17], [18], [33]. Conditional regularity results for 
weak solutions (of Serrin type) can be found in [9]. The smoothness of solutions arising 
from large axisymmetric data is addressed in [2] and [24J . Further regularity issues on the 
solutions have been discussed also [20[ [15] . 

The goal of this paper is to study in which way the variations of the temperature affect 
the asymptotic behavior of the velocity field. We point out that several different models are 
known in the literature under the name of "viscous (or dissipative) Boussinesq system". 
The asymptotic behaviour of viscous Boussinesq systems of different nature have been 
recently addressed, e.g., in [3J[T2]. But the results therein cannot be compared with ours. 

Only few works are devoted to the study of the large time behavior of solutions to (II. ip . 
See [21] [26]. These two papers deal with self-similarity issues and stability results for 
solutions in critical spaces (with respect to the scaling). On the other hand, we will be 
mainly concerned with instability results for the energy norm, or for other subcritical 
spaces, such as LP , with p < 3. 

A simple energy argument shows that weak solutions arising from data 8$ £ L 1 n L 2 
and uq £ l? a satisfy the estimates 



The above estimate for the temperature looks optimal, since the decay agrees with that 
of the heat kernel. On the other hand the optimality of the estimate for the velocity field 
is not so clear. 

For example, in the particular case 8$ = 0, the system boils down to the Navier- 
Stokes equations and in this simpler case one can improve the bound for the velocity into 
||^(i)||2 < 1 1 ito 1 1 2 - In fact, ||ii(t)||2 — > for large time by a result of Masuda [28]. Moreover, 
in the case of Navier-Stokes the decay of ||u(i)||2 agrees with the L 2 -decay of the solution 
of the heat equation. See [35], [25] [39] for a more precise statement. 

The goal of the paper will be to show that the estimate of weak solutions ||w(i)||2 < 
C(l + i) 1//4 can be improved if and only if the initial temperature has zero mean. To 
achieve this, we will establish the validity of the corresponding lower bounds for a class of 
strong solutions. 

In particular, this means that very nice data (say, data that are smooth, fast decaying 
and "small" in some strong norm) give rise to solutions that become large as t — > oo: our 
results imply the growth of the energy for strong solutions : 



u(t)\\ 2 <C(l + t) 



1/4 



and 



||^)|| 2 <C(l + t)- 3/4 - 



(1.2) 



c(l + t) 



/4 < ||^(i)||2 < CCl + t) 1 ^^ 



t » 1. 
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The validity of the lower bound in (jl.2p (namely, the condition c > 0), will be ensured 
whenever the initial temperature is sufficiently decaying but 

/ #o + 0. 

We feel that is important to point out here an erratum to the paper [23] . Unfortunately, 
the lower bound in (II, 2j) contredicts a result in [231 Theorem 2.3], where the authors 
claimed that [|u(£)[|2 — > under too general assumptions, weaker than those leading to our 
growth estimate. The proof of their theorem (in particular, of inequalities (5.6) and (5.8) 
in [23]) can be fixed by putting different conditions on the data, including J 9q = 0. This is 
essentially what we will do in part (b) of our Theorem 12 . 2 1 b elow . Similarity, the statement 
of Theorem 2.4 in [23] contredicts our lower bound (jl.6p below (inequalities (5.18)-(5.20) 
in their proof do not look correct). This will be also corrected by our Theorem 12.21 We 
would like to give credit to the paper [23] (despite the above mentioned errata), because 
we got from there inspiration for our results of Sections [3] and [H 

Our main tool for establishing the the lower bound will be the derivation of exact 
pointwise asymptotic profiles of solutions in the parabolic region \x\ ^> \ft. This will 
require a careful choice of several function spaces in order to obtain as much information 
as possible on the pointwise behavior of the velocity and the temperature. A similar 
method has been applied before by the first author in [6] in the case of the Navier-Stokes 
equations, although the relevant estimates were performed there in a different functional 
setting. 

Even though several other methods developped for Navier-Stokes could be effective for 
obtaining estimates from below (see, e.g., [13l [22j [30] ) , our analysis has the advantage of 
putting in evidence some features that are specific of the Boussinesq system: in particu- 
lar, the different behavior of the flow when \x^\ — > oo or when \J x\ + x\ — > oo, due to 
the verticality of the bouyancy forcing term Oe^ (see Theorem 12.61 below) . Moreover, the 
analysis of solutions in the region \x\ ;§> \ft and our use of weighted spaces completely 
explains the phenomenon of the energy growth: the variations of the temperature push 
the fluid particles in the far field; even though in any bounded region the fluid particles 
slow down as t — > oo (this effect is measured e.g. by the decay of the L°°-norm established 
in Proposition 12.50 . large portions of fluid globally carry an increasing energy during the 
evolution. Our result thus illustrates the physical limitations of the Boussinesq approxima- 
tion, at least for the study of heat convection inside fluids filling domains where Poincare's 
inequality is not available, such as the whole space. 

In fact, our method applies also to weighted L p -spaces, so let us introduce the weighted 
norm 

uu = (J \f(x)ni+\x\rd X 

Then we will show that strong solutions starting from suitably small and well decaying 
data satisfy, for t > large enough, 

(1.3) c(l + t)^ r+ l-V < \\u(t)\\ I *<C(l + tfr r+ l- 1 \ 
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for all 

3 

r > 0, 1 < p < oo, r H — < 3. 

P 

As before in (II. 2p . these lower bounds hold true with a constant c > as soon as the initial 
temperature has non-zero mean. Notice that in this case the L p -norms asymptotically blow 
up for large time if and only if p < 3. 

The above restriction r + | < 3 on the parameters is optimal. Indeed, under the 
same conditions yielding to (jl.3p we will also prove that when r>0, l<p<oo and 
t + | > 3, then one has [|u(i)||x,p = oo, for all t > 0. The fact that the L^-norm becomes 
infinite instanteneously for this range of the parameters is related to that the velocity field 
immediately spatially spreads out and cannot decay faster than |x| -3 as \x\ — > oo for t > 0, 
and this even if u £ Cfi°(R 3 ). 

As we observed before, the lower bound in (jl.3p brakes down when j 9q = 0. In such 
case, our decay estimates can be improved. We will establish the upper bound for weak 
solutions 

(1.4) lluWH^C'OL + t)- 1 / 4 

by the Fourier splitting method. This method was first introduced in [35] , We will need no 
smallness assumption on uq to prove (|1.4p . We have to put however a smallness condition 
of the form 

(1.5) ||0o||i<eo- 

We do not know if it is possible to get rid of fll .5 j) to establish (jl.4j) . Such smallness condi- 
tion however looks natural as it respects the natural scaling invariance of the system (jl.ip . 

As before, the decay estimate (|1.4p is optimal for generic solutions satisfying J 6q = 0. 
Indeed, when we start from localized and small velocity, we can establish the upper-lower 
bounds for strong solutions 

(1.6) c '(i + ^+!- 2 ) < \\ u{ t)\\ L? < c"(i + ^ (r+ §- 2) , 

for all 

3 

r > 0, 1 < p < oo, r H — < 4. 

p 

Similarily as before, the validity of the lower bound (the condition c' > 0) now requires a 
non vanishing condition in the first moments of 6. 

1.1. Notations. We denote by Cq° the space of smooth functions with compact support. 
The L p -will be denoted by || • |L and in the case p = 2 we will simply write || • ||. Moreover, 

V = {0GC o oo |V- ( /» = O}, 

denotes the completion of V under the norm || • and V the closure of V in Hq 
We will adopt the following convention for the Fourier transform of integrable functions: 

J-f(£) = /(£) = J f(x)e~~ 2nlx '^ dx. Here and throughout the paper all integral without 

integration limits are over the whole M 3 . 

The notations L p ' q and IX have a different meaning. For 1 < p < oo and 1 < q < oo, 

L p ' q denotes the classical Lorentz space. For the definition and the basic inequalities 
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concerning Lorentz spaces (namely the generalisation of the straightforward L p -L q Holder 
and Young convolution inequalities) the reader can refer to \27\ Chapter 2]. Here we just 
recall that L p,p agrees with the usual Lebesgue space L p , and L p '°° agrees with the weak 
Lebesgue space L v w (or Marcinkiewicz space) 

LP W = {/ : K" -> C, measurable , < oo}. 

The quasi-norm 

||/|| i& =supt[A / (t)]F 
t>o 

is equivalent to the natural norm on L p '°°, for 1 < p < oo. Here as it is usual we defined 

A/(s) = \{x : f(x) > s} 

where A denotes the Lebesgue measure. On the other hand, for 1 < p < oo and < r < oo, 
Lf is the weighted Lebesgue space, consisting of the functions / such that (1 + lxl) 9 / G L p . 
Notice that the bold subscript a introduced above is not a real parameter: in the notation 
Lct, this subscript simply stands for "solenoidal" . 

To denote general constants we use C which may change from line to line. In certain 
cases we will write C(a) to emphasize the constants dependence on a. 

We denote by e* A the heat semigroup. Thus, e tA UQ = J g t (x — y)uo(y)dy, where 
gt (x) = (4 7r )-3/2 e -|x| 2 /(4t) is the heat kernel. 

We denote by E(x) the fundamental solution of —A in M 3 . The partial derivatives of 
E are denoted by E x . , E x . Xk , etc. 

1.2. Organization of the paper. All the main results are stated without proof in Sec- 
tion [2j The statement of our theorems are splitted into two parts: part (a) is devoted to 
the properties of solutions in the general case (where the integral j 9q is not necessarily 
zero). In Part (b) of our theorems we are concerned with the special case J 6$ = 0. 

The rest of the paper is organized as follows. Sections [30 are devoted to the proof of 
our results on weak solutions and Sections [MS] to strong solutions. In Section [7] we collect 
a few technical remarks. 

2. Statement of the main results 

2.1. Results on weak solutions. We will consider weak solutions to the viscous Boussi- 
nesq equations and establish existence and natural decay estimates in LP spaces, 1 < p < 
oo, for the temperature, together with bounds for the growth of the velocities. Such esti- 
mates rely on the fact that in both equations of the system (jl.ip we have a diffusion term. 
They complete those of obtained in [16] where there was no temperature diffusion. 

We start recalling the basic existence result of weak solutions to the system See, 
e.g., PEG]. 

Proposition 2.1. Let (9o,uq) £ L? x L? a . There exists a weak solution (6,u) of the 
Boussinesq system il.l]) . continuous from M + to L? with the weak topology, with data 
(uq,8o) such that, for any T > 0, 

9 G L 2 (0, T; H 1 ) n L°°(0, T; L 2 a ), u G L 2 (0, T; V) n L°°(0, T; L 2 ). 
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Such solution satisfies, for all t G [0, T], the energy inequalities 

(2.1) ||^)|| 2 + 2 r||W(s)ll 2 ^<ll#ol| 2 . 

J o 

and 

(2.2) \\u{t)\\ 2 + 2 f\\Vu{s)\\ 2 ds < +c(||u || + * 2 |N| 2 
for all t > 0, and some absolute constant C > 0. 

One can improve the growth estimate on the velocity as soon as 6$ belongs to some LP 
space, with p < 2. For simplicity we will consider only the case 8o £ L 1 n L 2 . 

Moreover, it is natural to ask under which supplementary conditions on the initial 
data one can insure that the energy of the fluid ||u(t)|| 2 remains uniformly bounded. 
Theorem 12.21 provides an answer. Beside a smallness assumption on ||#o||i) we need to 
assume that f 6q = 0. 

In addition, it is possible to prove that not only remains uniformly bounded, but 

actually decays at infinity (without any rate). Explicit decay rates for ||it(i)|| can also be 
prescribed provided the linear part e A uo decays at the appropriate rate. 

Theorem 2.2. (a) Let (9q,uq) G L 2 x L 2 . Under the additional condition 9q G L 1 
the estimates on the weak solution constructed in Proposition \2.1\ can be improved 
into 



< C(t + l)-2 
(2-3) 

" i2 <c(t + r 



2 . 



Moreover, if 9q £ L 1 n LP, for some 1 < p < oo, then 

P <C(p)(t + l)-i (1 -^. 



(b) (The J 9q = case) In this part we additionally assume 8q G L\ and J 8q = 0. 
Then there exists an absolute constant eq > such that if 

(2-4) ||0 o ||i < eo 

then the weak solution of the Boussinesq system (jl.ip constructed in Proposi- 
tion \2.1\ satisfies, for some constant C > and all t G M + ; 

(2.5) ||0(t)|| 2 <C(l + t)-I 
and 

(2.6) ||u(t)|| 2 ^0 ast^oo. 

Moreover under the additional condition uq G L 3//2 n L 2 , we have 

(2.7) \\u{t)\\ 2 <c(i + ty 1/2 . 
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The proof of part (a) of Theorem 12,21 is straightforward. Part (b) is more subtle: its 
proof relies on an inductive argument: assuming that the approximate velocity -u ra_1 grows 
in L 2 at most like t 1//8 (which is actually better than provided by Proposition ^. ip we can 
prove that the same growth estimate remains valid for u n . A smallness condition on 9q is 
needed to insure that in the inequality ||u n (i)|| < C(l + i) 1 / 8 the constant can by taken 
independently on n. With this improved control on the growth of the velocity, applying 
the Fourier splitting method (introduced in |35j ) and a few boot-strapping we can improve 
our estimates up to the rates given in ()2.5[) and (|2.7D . 

Remark 2.3. Solutions of the Boussinesq system ([Lip with energy decaying faster than i -1 / 2 
might exist, but they are likely to be highly non-generic. Indeed, it seems difficult to 
construct such solutions assuming only that the data belong to suitable function spaces 
(possibly with small norms). The main obstruction is that one would need stringent can- 
cellations properties on the data, which however turn out to be non-invariant under the 
Boussinesq flow. A possible way to obtain such fast decaying solutions would be to start 
with data satisfying some special rotational symmetries, as those described in [5]. 

2.2. Results on strong solutions. The best way to prove the optimality of the estimates 
contained in Theorem 12.21 is to establish the corresponding lower bound estimates at least 
for a subclass of solutions. For the study of the estimates from below we will limit our 
considerations to a class of strong solutions. This is not a real restriction as lower bound 
estimates established for solutions emanating from well localized, smooth and small data 
are expected to remain valid in the larger class of weak solutions. Studying strong solutions 
has also the advantage of better putting in evidence some interesting properties specific of 
the Boussinesq system, such as the influence of the vertical buoyancy force on the pointwise 
behavior of the fluid in the far-field. 

The existence of strong solutions to the system (II. ID will be insured by a fixed point 
theorem in function spaces invariant under the natural scaling of the equation. Thus, if 
u G X , where X is a Banach space to be determined, we want to have, for all A > 0, 

II^aIIa' = IMIaTj where u\(x,t) = Xu(Xx, X 2 t) 

is the rescaled velocity. A suitable choice for norm of the space X, inspired by |10j . is 

(2.8) \\ u \\x = ess sup (Vt + \x\)\u{x,t)\. 

xem. 3 ,t>o 

This choice for X is quite natural. Indeed, whenever |«o(a:)| < C|x| -1 , the linear evolution 
e tA uo belongs to X and this paves the way for the application of the fixed point theorem 
in such space. 

More precisely, we define X as the Banach space of all locally integrable divergence-free 
vector fields u such that \\u\\x < oo, and continuous with respect to t in the following 
usual sense: u(t) — > u(0) in the distributional sense as t — > and ess sup^g^ \x\ \u(x,t) — 
u(x,t')\ ->• as t t' if t' > 0. 

In the same way, if 9 belongs to a Banach space y, we want to have 

||0 A ||3; = \\e\\ y , where 9\(x, t) = \ 3 6(\x, X 2 t) 



s 
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is the rescaled temperature. We then define a Banach space y of scalar functions through 
the norm 

(2.9) ||0|| y = ||6»|| L oo (L i) +ess sup (y/t + \x\f\9(x, t) \ 

' ieR 3 ,t>o 

and the natural continuity condition on the time variable as before. 

The starting point of our analysis will be the following proposition providing a simple 
construction of mild solutions (u, 9) € Xxy. We will refer to them also as strong solutions. 
Indeed, one could prove that such solutions turn out to be smooth, as one could check by 
adapting to the system classical regularity criteria for the Navier-Stokes equations 
like that of Serrin |37j . See the paper [9]. The smoothness of these solutions, however, 
plays no special role in our arguments. 

Proposition 2.4. There exists an absolute constant e > such that if 

(2.10) ll#o||i < e> ess sup |x| 3 |0o(^)| < e> ess sup |x| |txo(^)| < e 

where uq is a divergence- free vector field, then there is a constant C > and a (mild) 
solution (u,9) £ X x y of (jl.ip . such that 

(2.11) HuH* < Ce and \\6\\y < Ce. 
Moreover, these conditions define u and 6 uniquely. 

Next Proposition shows it is possible to obtain better space-time decay estimates pro- 
vided one starts with suitably decaying data. 



Proposition 2.5. (a) Let uq and 6$ as in Proposition \2.4\ and satisfying the addi- 
tional decay estimates, for some 1 < a < 3, b > 3, and a constant C > 0, 

\u (x)\ < C(l + |x|)" a , 
(2-12) , 

i0oOzoi<c(i+i*ir 6 . 

Then the solution constructed in Proposition \2.J\ satisfies, for another constant 
C > independent on x and t, 

(2.13) \u(x,t)\<C inf + t)^- 1)/2 

0<ri<a 

and 

(2.14) \0{x,t)\<C inf \x\~^{l + t)^' 2 . 

0<r)<b 

(b) (the j 9q = case) Assume now 2 < a < 4, a / 3, and b > 4 and let uq and 9q 
satisfying the previous assumptions. If, in addition, 

(2.15) J 9 = and 9 G L\, 

then the decay of u and 9 is improved as follows: 

\u(x,t)\<C inf (1 + |x|)- r? (l + i)^- 2 >/ 2 , 

0<r;<a 

(2.16 " , „ 

\9(x,t)\<C ini ■ (l + |x|)-"(l + t)^)/i. 

0<ri<b 



DECAY AND GROWTH OF SOLUTIONS OF BOUSSINESQ (July 29, 2010) 9 

Recall that the fundamental solution of —A in IR 3 is E{x) = c|x| _1 . Thus, in the 
following asymptotic expansions, VE X3 and VE XhtX3 are vectors whose components are 
homogeneous functions of degree —3 and —4 respectively. In particular \VE x . 3 (x)\ < 
C\x\~ 3 and \VE x . ;X3 (x)\ < C|x|" 4 . 

We are now in the position of stating our main results on strong solutions. The first 
theorem describes the asymptotic profiles of solutions in the parabolic region |x| ;§> \ft. 
Roughly, it states that all sufficiently decaying solutions (u, 9) of (jl.ip behave in such 
region like a potential flow. 

Theorem 2.6. 

(a) Let a > | and b > 3. Let (u,9) be a (mild) solution of (jl.ip satisfying the decay 
estimates (|2. 13|) - (|2. 14|) . Then the following profile for u holds: 

(2.17) u(x,t) = e tA u (x) e^jt (VE X3 )(x) + K{x,t) 

where 1Z(x,t) is a lower order term with respect to tX7E X3 (x) for \x\ ;§> \ft, 
namely, 

(2.18) lim ^4 = 0. 

Vt 

(b) (the J 9q = case) Assume now a > 2 and b > 4. Assume also that J 9q = 0. 
Let (u, 9) be a solution satisfying the decay condition (|2.16p . Then the following 
profiles for Uj (j = 1, 2, 3) hold: 

(2.19) Uj (x,t) = e tA u {x) -VE XjX3 (x) ■ ( [ [ y9(y,s)dyds) +K(x,t) 



where 1Z is a lower order term for \x\ S> ^ft >■ 1, namely 

(2.20) lim ^4 = 0. 

f,M^ t\x\~* 
Vt 

The following remark should give a better understanding of the theorem. 

Remark 2.7. (a) (The case j 9q ^ 0). We deduce from the asymptotic profile (I2.17|) 
the following: when \e tA uo(x)\ <C i|x| -3 (this happens, e.g., when we assume also 
\uo(x)\ < C|x|~ 3 and \x\ > \ft » 1) and / 9 ^ then 

(2.21) u(x,t) ~ (^J 9^j t (VE X3 )(x), for \x\ > v 7 *- 

(The exact meaning of our notation and of statements (|2.2ip and (|2.22p below is 
made precise in the proof), 
(b) (The / 9q = case) We deduce from the profile (I2.19P the following: when 
\e tA uo(x)\ <C i|x| -4 (this happens, e.g., when we assume also |uo(x)| < C|x| -4 
and |x| » \/t 3> 1) then 

(2.22) Uj (x,t) ~ -VE XjX3 (x) ■ (J J y0(y,s)dyds\ for \x\ » v 7 * > 1. 
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A remarkable consequence of the previous theorem is the following. 
Corollary 2.8. 

(a) Let a > |, b > 3 and let (u,9) be a solution as in Part (a) of Theorem \2.b\ Then 
for all r, p such that 

3 . . 

T > 0, 1 < p < oo, r H — < minja, 3j, 

p 

there exists to > such that the solution satisfies the upper and lower estimates in 
the weighted-LP -norm 

(2.23) 0(|m o |)(l + t)* (r+ S- 1) < Ht)|| i? <C' / (l + t)' (r+ *~ 1) 

for all t > to- Here, mo = J 9o and <p: R + — > R is some continuous function such 
that (f)(0) = and <f>(a) >0 if a >0. 

(b) (the J 9o = case) Under the assumptions of the previous item, with the stronger 
conditions a > 2, b > 4 and the additional zero mean condition mo = 0, let us set 
m = liminf^oo j f^f y9(y, s) dyds . Then, for all r,p such that 

3 

r > 0, 1 < p < oo, r H — < minja, 4}, 

p 

we have 

(2.24) (f)(m) (l + i)5 (r+ t- 2) < ||n(t)|| L P < C'(l + t) 5 (r+ l" 2) 

for another suitable continuous function (f>: M + — > R suc/i t/iat (f>(0) = and 
^>(a) > for a > 0. 



Remark 2.9. When J" # 7^ 0, we thus get by (|2.23p the sharp large time behavior ||u(i)|| L p ~ 

When J 9o = and m/Owe have the faster sharp decay — t^ r+ p 2 . The 

condition m / is satisfied for generic solutions. It prevents 9 to have oscillations at 
large times. 

3. The mollified Boussinesq system and existence of weak solutions 

The existence of weak solutions to the Boussinesq system is well known, see [9]. Their 
uniqueness, however, is an open problem. Moreover, we do not know if any weak solutions 
satisfy the energy inequality and the decay estimates stated in Proposition 12.11 For this 
reason, we now briefly outline another construction of weak solutions, which is well suited 
for obtaining all our estimates. 

We begin by introducing a mollified Boussinesq system. As the construction below is 
a straightworward adaptation of that of Caffarelli, Kohn and Nirenberg, [8], we will be 
rather sketchy. For completeness we recall the definition of the "retarded mollifier" as 
given in [8]. Let ip(x,t) G C°° such that 



4>> 0, 



f'OQ f' 

J Jipdxdt = l, suppV" C {(x,t) : \x\ 2 < t, 1 < t < 2}. 



DECAY AND GROWTH OF SOLUTIONS OF BOUSSINESQ (July 29, 2010) 



11 



For T > and u G L 2 (0,T; L\ ), let u : R 3 x R ->■ R 3 be 

'u(ar,t) if (x,t) G R 3 x (0,T), 
otherwise. 



Let (5 = T/n. We set 

^ 5 (n)(x, i) = 6~ A J "^{jji -g) u(x -y,t-r) dydr. 

Consider, for n = 1,2,... and 6 = T/n, the mollified Cauchy problem 
' d t 6 n + ^s(u n ~ l ) ■ V9 n = A9 n 

(3.1) I d t u n + V • (^ 5 (u n ) u n ) + Vp n = Au n + 9 n e 3 x G R 3 , t G R+ 

, V • -u n = 0. 

with data 

(3.2) 9 n \ t =o = e and u n | 4=0 = u . 

The iteration scheme starts with u° = 0. 

Note that since div-u = 0, we also have div (^s(u n )) = 0, for t G R+. At each step n, one 
solves recursively n + 1 linear equations: first one solves the transport-diffusion equation 
(with smooth convective velocity) for the temperature; after 9 n is computed, solving the 
second of (|3.ip amounts to solving a linear equation on each strip R 3 x (m6, (m + 1)6), 
for m = 0, 1, . . . , n — 1. 

For solutions to (13. ip we have the following existence and uniqueness result, 

Proposition 3.1. Let (9q,uq) G L 2 x L 2 . For eac/i n G {1,2,...} £/iere exists a unique 
weak solution (9 n ,u n ,p n ) of the approximating equations with data \3. Sty such that, for 
any T > 0, 

6 n G L 2 (0, T; H 1 ) n L°°(0, T; L 2 ), u n G L 2 (0, T; V) n L°°(0, T; L 2 ) 

and 

p n G l 5 / 3 (0,T;L 5 / 3 ) + L°°(0,r;L 6 ) 

Moreover, for all t > 0, 6* n and u n satisfy the energy inequalities as in (12. ip and (I2.2p . 
in particular, the sequences 9 n , u n and p n , n = 1,2,... are bounded in their respective 
spaces. 

Proof. 

This can be proved using the Faedo-Galerkin method. As the the argument is standard 
(see, e.g. [38l Theorem 1.1, Chapter III], or [8j Appendix]), we skip the details. We only 
prove the condition on the pressure since this is the only change that we have to make 
in®. 

Taking the divergence of the second equation in (13. If) we get p = p™ + p%, where 

Ap^ = -J2did 3 Ku]) 

and 

Ap n 2 = d X3 9 n . 
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Thus, p" G -L 5 / 3 (0, T; L 5 / 3 ) uniformly with respect to n, by the energy inequality for u n , 
interpolation, and the C alder on- Zygmund theorem, as proved in [8j. On the other hand, 
— P2 = E Xa * 9, where E(x) is the fundamental solution of —A. Thus, E X3 (x) = 

belongs to the Lorentz space 

L 3/2,oo( R 3)_ But gn £ L°°(0,T;L 2 ) uniformly with respect 
to n, hence Young convolution inequality in Lorentz spaces (see [27J Chapter 2]) yields 
P2 G L°°(0, T;L 6 ) uniformly with respect to 11. 

□ 

It follows from Proposition 13.11 that, extracting suitable subsequences, p n = p\ + p%, 
where (p™) converges weakly in L 5 / 3 (0, T; L 5 / 3 ) and P2 converges in L°°(0, T;L 6 ) in the 
weak-* topology. Moreover, (u n ) is convergent with respect to the topologies listed in [HI p. 
828-829]. On the other hand (8 n ) will be convergent with respect to the same topologies, 
because all estimates available for u n hold also for 9 n . 

No additional difficulty in the passage to the limit in the nonliner terms arises, in the 
equation of the temperature other those already existing for the Navier-Stokes equations. 
Hence, the distributional limit (8,u,p) of a convergent subsequence of (9 n ,u n ,p n ) is a 
weak solutions of the Boussinesq system. This establishes Proposition 12.11 

□ 

We finish this section by establishing the natural L p -estimates for the approximating 
temperatures: 

Lemma 3.2. Let (9q,uq) G L 2 x l?^ and let (9 n ,u n ,p n ) be the solution of the mollified 
Boussinesq system (|3.ip for some n G {1, 2, . . .}. Let also 1 < p < 00. If 9o £ L 1 Pi L p , 
then 

ct . A-K 1 -?) 



(3.3) l|0 n (*)llp< \\0o\\i{j + A 



where A = A(p, ||#o||i) ll^ollp) an d c > is an absolute constant. 

Proof. First notice that, for each n, 9 n is the solution of a linear transport-diffusion 
equation with smooth and divergence- free velocity ^(u n ~ 1 ). The LP decay estimates for 
these equations are well known. We reproduce the same proof as in |14| 119] ) expliciting 
better the constants, as we will need the expressions of such constants later on. 
A basic estimate (valid for 1 < p < 00) is 

(3-4) P n (t)\\ P < 

See [14, Corollary 2.6] for a nice proof of (j3.4|) that remains valid in the much more general 
case of trasport equations with (or without) fractional diffusion. 

We start with the case 2 < p < 00. Multiplying the equation for 9 n by p\9 n \ p ~ 2 9 n and 
integrating we get 

±\\9^t)\\i + iiv(rr/ 2 )(t)n 2 < 0. 

at F p 

By the Sobolev embedding theorem, H 1 C L 6 , hence 

\\e n (t)\\ p 3p <c\\v(\9T /2 )(t)\\ 2 - 
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The interpolation inequality yields 

n# n (t)ii P < iie ft (t)ii; /( *" l) iifl ri (t)iig , " l)/( *" 1) . 

Combining these two inequalities with the basic estimate ||# n (i)||i < ll^olli) we obtain the 
differential inequality 



dt 

Integrating this we get 



i r(t)llg -" c P ii6biiy /( *- 3) (r(t)llg) 



n (oiif;< 7^777^7 + 



t i \ -3(p-l)/2 



p ~ \3C P \\e \\f p)/ ^- 3) |NI? p)/(3p ~ 3) 

and estimate (|3.3|) follows with 

(2p)/(3p-3) 



/I 



7 0||1 



The case 1 < p < 2 is deduced by interpolation. 
In the p = 2 case, we obtain the following 



□ 



Lemma 3.3. Let 9q G L 1 nL 2 anduo G L 2 . Lei (9 n ,u n ,p n ) be the solution of the mollified 
Boussinesq system (13. ID /or some n G {1, 2, . . .}. Then, 

\\e n (t)\\<\\Oo\\i(ct + A )- 3/4 , 
\\u n (t)\\ < WuoW + c'weoht 1 ^, 

for two absolute constants C, C > and j4o = (H^olli/H^olh) 4 ^ 3 - 

Proof. We only have to estimate the L 2 norm of the velocity. We make use of the identity 

^K(*)lli = / 

that can be justified exactly as for the mollified Navier-Stokes equations, see [8| and [55] . 

Multiplying the velocity equation in the Boussinesq system (|3.ip by u n and integrating 
we get 

(3.6) 4lK(*)ll 2 + 2||Vn"(t)|| 2 < 2\\u n 

dt 



Dividing by \\u r 



ih n (t)\\ < 



dt' 

Now we use the decay of ||# n (t)|| obtained in Lemma [3. 2 i Integrating we obtain the second 
of (HLH). 

□ 
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4. Improved bounds for weak solution in the case J 9 = 

The estimates obtained in Lemma 13.31 can be considerably improved provided we ad- 
ditionally assume J 9q = and the moment condition 9q G L\. First of all, from an 
elementary heat kernel estimate one easily checks that in this case 

(4.1) ||e' A 0o|| 2 <-4i(i + ir 5/2 , 

where Ax > depends only on the data, through its L 2 -norm and the J \x\ \9o(x)\dx 
integral. 

We will now see that in this case the approximated temperature 9 n (t) also decays at 
the faster rate (t + 1) _5//4 in the L 2 -norm. Using this new decay rate for 9 n it is possible 
to show that the velocities are uniformly bounded in L 2 . Once we have a solution such 
that ||« n (t)|| remains bounded as t — > oo one can go further and prove that u n actually 
decays at infinity in L 2 at some algebraic decay rate, depending on the decay of the linear 
evolution e <A uo- In view of the passage to the limit from the mollified system (|3,ip to the 
Boussinesq system (jl.ip . all the estimates must be independent on n. 

Proposition 4.1. Let (9q,uq) £ (L\ n L 2 ) x L 2 and assume J 9q = 0. There exists an 
absolute constant £o > such that if 

(4.2) ||0 o ||i < eo 

then the solution of the mollified Boussinesq system \3. with data (uq,9q) satisfies 

(4.3) ||^(t)|| 2 <A(l + t)~f 
and 

(4.4) \\u n (t)\\ 2 <A 

for all n £ N and t G M + . Here A > is some constant depending on the data uq and 9q, 
and independent on n, and t. 

Proof. We denote by C a positive absolute constant, which may change from line to 
line. We also denote by Ax, A^, . . . positive constants that depend only on the data. More 
precisely, Aj = Aj (||0 O ||, ||0o||li, IKH) • 

Step 1: An auxiliary estimate. 

We make use of the Fourier splitting method introduced in [35 j . The first step consists 
in multiplying the temperature equation by 9 n and to integrate by parts. Using the 
Plancherel theorem in the energy inequality for 9 n , we get 

\jt j l^)! 2 ^- / lel 2 !^)! 2 ^ 
Now split the integral on the right hand side into 5U5 C , where 
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and A; is a constant to be determined below. Noting that for £ E S c one has — |£| 2 
— 2 (t+i) » ^ follows that 



d 
dt 



\e n (t,t)\ 2 < 



t + i 

k 

t + i 



\e^,t)\ 2 

c 

r(£,t)i 2 d£+ 



t+i 



\e n &t)\ 2 <%. 



Multiplying by (1 + t) k we obtain 



(4.6) 



d_ 
~dt 



/ r(^)i 2 d£ 



^fc^+i)*- 1 / i^(e,t)i 2 de, 



where S is as in (|4.5p . 

Taking the Fourier transform in the equation for # n in (|3,ip we get 



(4.7) 
Hence 



|^,t)| 2 <2| e "^ 2 ^ | 2 + 2|£| 



V^tollirtoll* 



Replacing this in (|4.6p and applying the Plancherel theorem we get 



d 

dt L 



(1 + tfF 



feimn/i\||2 



< c 



„/An i|2 



[|2 (1+t jfc-l + (1+t jfc-7/2 / \\ u ^\\\\0n\\ ds 



From where it follows, letting k = 7/2, 



\\9 n (t)\\ 2 <(l + tr 7 / 2 \\\9of + C 



(4.8) 



f\\e^e \\ 2 {l + sf' 2 ds 



+ 



\ 2 

ll^WIIII^WUdr ds 



10 vo 

Recalling the estimate (|4.ip for the linear evolution we get, for all n E N, 



(4.9) ||0 



n/ + \||2 



< c 



Al( l + r 5/2 +(1 + r 7/ 2 £ ^ 



u n ~ (r)\\ \\9 n (r)\\dr) ds 



We now use the following inequality, deduced from estimate (|3.5p . 

(4.10) \\9 n (t)\\ <C7||# ||i^ 3/4 . 

Putting this inside (|4.9p we obtain a new bound for ||# ra || 2 , namely 
(4.11) 



\\9 n (tW<C 



A 1 (l + t)- 5 / 2 + \\9 \\l(l+t)- 7 / 2 J* QT 



IvT' 1 (r) II r~ 3 / 4 dr ds. 



Step 2: The inductive argument. 

We will now prove by induction that, for all positive integer n we have 
(4.12) IK^WII < WuoW + Mt 1 / 8 
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where M > is some constant independent on n (but possibly dependent on the data 9q, 
uq) to be determined. Notice that estimate (I4.12p is actually better than what we have so 
far (compare with the second of (|3.5p ). 

For n = 1 the inductive condition (|4.12p is immediate since u° = 0. Let us now prove 
that ||tt n || < \\u \\ + Mt 1 / 8 assuming that (|4,12p holds true. 

We get from (|4.1ip and the induction assumption (|4.12p 

12 < C A x (l + 1)~ 5 / 2 + [|^o[|?||i*o|| 2 (l + 1) -2 + Af 2 |[» ||?(l + 1) _7/4 



This implies 
(4.13) 



< C 



V^i + INIiflKII + M 



-7/8 



This last inequality will be used to estimate ||i£ n || as follows. First recall that 

d 



(4.14) 



dt 



\u n (t)\\ < \\e n 



After an integration in time we get, using (|4.13j) . 



\u n (t)\\ < iKH + c 



v^4i + ||^o||if||«oll + M 



,1/8 



For the induction argument we need to prove ||ti n || < ||ito|| + 

Mt 1/8 . Hence we need that 



C 



v^i + ||0o||i(lKII +m) 



< M. 



Choosing M large enough, for example, 

M = max{||ii ||,2CV r Ai}, 
our condition then boils down to the inequality 

||0o||i < 1/(4C). 

The validity this last inequality is insured by assumption (|4.2p . This concludes the induc- 
tion argument and establishes the validity of the estimate (|4,12p for all n. 

Step 3: Uniform bound for the L 2 -norm of the velocities u n . 

The result of Step 2 implies the existence of a constant A2 > such that, for all n > 1, 



(4.15) IK -1 (i)H 2 <^ 2 (l + t) 1/4 . 

In the proof of the previous Step we also deduced that, for some ^3 > 0, 

l|r(t)|| 2 <^ 3 (i + t)- 7/4 . 

Combining such two estimates with inequality (|4.9p we easily get 

\\0 n (t)f <Ai(l + t)- 2 . 
Now using this improved estimate for ||# n || 2 with (|4.15p in (|4.9|) arrive at 

\\9 n (t)\\ 2 < A 5 (l + t)- g / A . 
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Going back to the differential inequality (|4.14p we finally get, for some constant Aq > 
independent on n, and t £ R + , 

\\u n (t)\\ 2 <A 6 . 

Replacing in (I4.9P we can further improve the decay of 9 n up to 

\\e n (t)\\ 2 <A 7 (i + ty 5 / 2 . 

□ 

We now would like to improve the result of the previous Proposition by establishing 
decay properties for ||u ra (i)|| 2 . Specifically, if we assume in addition that the linear part 
of the velocity satisfies ||e* A iio|| 2 < C(l + i) -1 / 2 (this happens e.g. when u G L 3 / 2 n L 2 ), 
then the same decay holds for the approximate velocities u n . 

Proposition 4.2. Let (9o,uq) € (L\nL 2 ) x L 2 ). Assume also that J 9$ = and ||#o||i < 
Eq, where Eq is the constant obtained in the previous proposition. 
Then the approximate solutions of (|3.ip satisfy, 

||-u n (t)|| ->• 0, as t -> oo, 

uniformly with respect to n. Moreover, if uq E L 3//2 n L 2 a , then 

(4.16) \\u n (t)\\ 2 < Ail + t)- 1 / 2 , 
for some constant A > independent on n, and t. 

Proof. We denote by A > a constant depending only on the data that might change from 
line to line. The proof follows by Fourier Splitting. Since the estimates are independent 
of n, we simply denote the solutions by (6,u). Multiply the second equation in (13. ID by 
u, integrate in space. By Proposition 14.11 we get 

(4.17) 4lK*)H 2 + 2||Vn|| 2 < C(t + l)~f 

at 

Arguing as for the proof of inequality (|4.6p we obtain 

< Ck{t + l) fc_1 J \u{i, t)\ 2 d£ + C(t + l) fc "t 

Where S was defined in (I4.5|) . From now on, k = 7/2. 

We need to estimate |u(£, t)\ for £ £ S. Computing the Fourier transform in the equation 
for u in (|3.ip . next applying the estimate ||n(t)|| 2 < A obtained in Proposition 14. II we get 

m,t)\ <e-*l«l 2 |n5| + |e| f \\u(s)\\ 2 ds + fie^s^ds 

Jo Jo 
ft 

<e- t ^ 2 \-u^\+At\C\+ / \0(£,s)\ds. 

Jo 

But computing the Fourier transform in the equation for 9 in (|3.ip and applying once 
more the estimates of Proposition 14.11 we have 

\m,s)\<e- s W 2 \9 \+A\Z\<\9 (O\+A\Z\ 

Hence, 

(4.19) \u(tt)\ 2 <A 



(4.18) 



d 
~dt 



(t + i) fc K*)[| s 



e- 2 ^ |2 (N 2 + i^oi 2 )+t 2 iei 2 ■ 
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Integrating on S and applying inequality ()4. 1[) we deduce 

| e ^ f + (l + r V2 



J s \u(i,t)\ 2 dt<A 



Putting this inside (14.180 and integrating on an interval of the form [t e ,t], with e > 
arbitrary and t e chosen in a such way ||e tA uo|| 2 < e f° r t > t e , we obtain 

||«(t)|| — > 0, as i — >• oo. 

On the other hand, in the case no £ L 3 / 2 n L 2 , we have 

^|n(e,t)| 2 ^<^(l + t)- 1/2 . 

Now going back to (|4,18p and integrating in time we finally get 

\\u(t)\\ 2 < Ail + ty 1 ' 2 . 

□ 

We are now in the position of deducing our result on weak solutions to the Boussinesq 
system (jl.ip . 



Proof of Theorem 12.21 Now this is immediate: passing to a subsequence, the approxi- 
mate solutions 9 n and u n converge in L 2 oc (R + , M 3 ) to a weak solution (9, u) of the Boussi- 
nesq system Moreover, the previous Lemmata imply that 9 n and u n satisfy estimates 
of the form 



l« B («)ll </(*). 



for all t > 0, 



where f(t) is a continuous function independent on n. Then the same estimate must 
hold for the limit 9 and u, except possibly points in a set of measure zero. But since weak 
solutions are necessarily continuous from [0, oo) to L 2 under the weak topology, ||0(i)|| and 
!!«(£) || are lower semi-continuous and hence they satisfy the above estimate for all t > 0. 
This observation on the weak semi-continuity is borrowed from |25j . 

□ 



5. Strong solutions: preliminary lemmata 
The integral formulation for the Boussinesq system, formally equivalent to (jl.lB reads 



(5.1) 



0(t) = e tA 9 - / e (*-*) A V • (9u)(s) ds 
Jo 

u(t) = e tA u - [ e ('- s ) A PV • (u®u)(s)ds+ [ 
Jo Jo 



P9(s)e 3 ds. 



{ V • u = 



The above system will be solved applying the following abstract lemma, which slightly 
generalizes that of G. Karch et N. Prioux (see [26\ Lemma 2.1]). 
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Lemma 5.1. Let X and y be two Banach spaces, let B: X x X — >• X and B : y x X — >• y 
be two bilinear maps and L: y — )• X a linear map satisfying the estimates \\B(u,v)\\x < 
a!i[|^[|^'||v||Afj \\B(0, v)\\y < oc2\\0\\y\\u\\x and \\L(9)\\x < a 3||$ll;y; f or some positive con- 
stants ax, a.2 and a^. 

Let < rj < 1 be arbitrary. For every (U, 0) G X x y such that 

t? P7 \\x +«3 6 y < 77. : 

(2«i + a 2 ) 

i/ie system 

(5.2) = 9 + 5(0, it), it = C/ + B(u,u) +L(0) 

/ias a solution (u,9) £ X xy. This is the unique solution satisfying the condition 

r]\\u\\x + oi3\\9\\y ^ V(l - v)( 2a i + a 2)- 

Proof. In the case < < 1, one can take r] = 03. In such particular case, this lemma 
is already known, see |26t Lemma 2.1]. Therefore, we only have to prove that we can get 
rid of the restriction < 03 < 1. This is straightforward. We introduce on the space y 
an equivalent norm, defined by \\0\\y = 7fr\\9\\y- By Lemma 2.1 of Karch and Prioux, 
applied in the space (X, || • \\x) and (y, \\ • \\y), we have that if 

\\ U \\X + O y ^ 77, 1 y2~, 

then the system f)5.2[) has a unique solution such that 

\\ u \\x + \\Q\\y < (1 - f?)(2«i + a 2 ). 
The conclusion of Lemma 15. II is now immediate. 

□ 



Remark 5.2. Using our improved version of Lemma 2.1 of [26], it is be possible to get rid 
of the smallness assumption |/3| < 1 in the main results of Karch and Prioux [26] . 

Remark 5.3. The proof of Lemma 2.1 in |26j relies on the contraction mapping theorem. In 
particular, the solution can be obtained passing to the limit with respect to the X x^-norm 
in the iteration scheme (k = 1, 2, . . .) 

(u°,9°) = (U,@), 

(5 3) 

(u k+1 ,9 k+1 ) = (u° + B(u k ,u k ) + L(9 k ),9° + B(u k ,9 k )). 
See [26] for more details. See also [331 Lemma 4.3] for similar abstract lemmata. 

Let a > 1. We define X a as the Banach space of divergence-free vector fields u = u(x, t), 
defined and measurable on R 3 x IR + , such that, for some C > 0, 

(5.4) \u(x,t)\<C inf \x\- v {l + t) (r '- 1)/2 . 

0<ri<a 

In the same way, for b > 3 we define the space 3^ of functons 9 E L^°(L 1 ) satisfying the 
estimates 

(5.5) \9(x,t)\<c inf bnn+t)^- 3 )/ 2 . 

0<Ty<6 
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Such spaces are equipped with their natural norms. 

They are obviously decreasing with respect to inclusion as a and b grow. Recalling the 
definition of X and y in Section 12.11 we see that X\ = X n L^ t with equivalence of the 
norms and that 3^ = y H L^ t . 

The estimates (|2.13|) - ()2.14j) in Proposition 12.51 are thus equivalent to the conditions u G 
X a and 9 £ 3V 

We start with some elementary embeddings. 

Lemma 5.4. Let L p,q be the Lorentz space, with 1 < p < oo and 1 < q < oo. Then the 
following four inequalities hold: 



1/3 



LP'i 



< C\\u\\ x t* { p 1] , 3 < p < oo, 



, . Mt)\\ LP>q <C\\u\\ Xa (l + t)^l 1) , f<p<oo, 

(5.6) j 3 

||^(t)||L^ < C\\9\\y t^p- 3) , 1< p < OO 

||0(t)[|zP,« < C||0b 3 (l + t) |( ^~ 3) , 1 < p < oo 
/or some constants C depending only on p and q. In particular, choosing p = q one gets 
the correponding estimates for the classical L? -spaces. 

Proof. The above estimates for the weak-Lebesgue spaces L p '°° are simple. Indeed, if 
u 6 X, then |u(x,t)| < C|x| -3 / p t2(p -1 ) anc j one nas only to recall that any function 
bounded by |x| _3//p belongs to L p '°°. The other L p '°°-estimates for u and 9 contained 
in ()5.6p In the case 1 < q < oo, we use that L p ' q is a real interpolation space between 
£P-e,oo an( j jj>+£,co ^ £ Qr gji i <p_£<p<p-j_ e <oo. Therefore estimates (|5.6p for all 
1 < q < oo follow from the corresponding estimates in the particular case q = oo via the 
interpolation inequality. 

□ 

The first useful estimate in view of the application of Lemma 15.11 is the following. 

Lemma 5.5. Let 1 < a < 3 and 9 £ 3^3 • We have, for some constant C > depending 
only on a, 

(5-7) \\L(9)U a < C\\9\\y s . 

Moreover, 

(5.8) ||L(0)[|*<C||0b. 



Proof. We prove only (|5.7p since the proof of (]5.8p is essentially the same. By a renor- 
malization, we can and do assume that \\9\\y 3 = 1. Let M.(x, t) be the kernel of the 
operator e* A P. Then we can write 



L(9)(x, t) = J jK(x-y,t- s)9(y, s)e 3 dy ds. 

We have the well known estimates for K (see, e.g., Prop. 1]) 
(5.9) |K(x, t)\ < C\x\- a r { ?- a W 2 , for all < a < 3 
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where C > is come constant independent on x,t and < a < 3. We also recall the 
scaling relation 

(5.10) K(x, t) = t-^ 2 K(x/Vi, 1) 

and the fact that K(-, t) € C°°(M 3 ) for t > 0. The usual L p estimates for K are 



(5.11) \\K(t)\\ p < Ct~ 3/2+ ^ , Kp<oo. 



2 r2 



Using the L -L convolution inequality, we get 

||£(0)(t)lloo < C At - s)- 3 /\l + s)-VUs < C(l + i)~ 1/2 - 



Owing to this estimate, the conclusion -L(#) G <f a will follow provided we prove the point- 
wise inequality, 

\L(6)\(x,t) < C a \x\- a t^~^l 2 , V(x,t) s.t. |x|>2 v / t. 
This leads us to decompose 

L(0) = /i+/ 2 + / 3 , 

where Ji = JoXl<M/2 • • • ' h = Sof\x- v \<\x\/2 ~- and / 3 = JoI\y\>\ x \/ 2> \x-y\>\x\/2 • • • ■ Usin S 

€ Lf> (L 1 ) we get 

(5.12) |Ji|(M) < C|x|- 3 t, 

which is even better in the region {(x,t): \x\ > 2\ft\ than what we need (recall that 

1 < a < 3). Using now |#(x,i)| < C| as| — 3 and the scaling properties of K we obtain by a 
change of variables 



\h\(x,t) <C\x\' 3 [ [ 

JoJ\y 

< C\x\~ 3 t \og{\x\ /yft) 



\K{y,l)\dy 

\<\x\/{2j!) 



< c 1 I i-o^fq-iya 

— 3— a I I 

for |x| > 2y/i, and 1 < a < 3. Next, using again \6(x,t)\ < C|a;|~ 3 and |K(sc,i)| < C|x|~ 3 , 

\h\(x,t)<c[ [ \y\' 6 dy ds < C\x\~ 3 t. 

JaJ\ y \>\x\/2 

Therefore, 

\L(6)\(x,t) < ^x^H^-VI 2 , \x\ > 2y/i. 
Lemma 15.51 in now established. 

□ 



We collect in the following Lemma all the estimates on B(u, v) that we shall need. (We 
will apply estimate (|5.15p in the proof of Proposition 12.41 estimate ()5.14j) for Proposi- 
tion (|2.4p and estimate (|5.14p in Theorem 12.61 



22 



LORENZO BRANDOLESE AND MARIA E. SCHONBEK 



Lemma 5.6. Let 1 < a < 3. For some constant C > 0, depending only on a we have 

(5.13) ||-S(u> v)\\x — ^\\ u \\x\\ v \\x 
and 

(5.14) \\B(u,v)\\ < C\\u\\ \\v\\ , 

where (2a)* =min{2a, 4}. Moreover, 

(5.15) ||S(u,v)||y < C • 

Proof. We begin with the proof of the first estimate. As before, we can assume \\u\\ x = 
U L =1. We start writing 

II iiri a 

(5.16) B{u,v)(x,t) = j j F(x — y,t — s)(u ® v)(y, s) dy ds, 

where F(x,t) is the kernel of the operator e tA FV . 

The well known counterpart of relations (|5.9p - (|5.10p are (see, e.g., [29], Prop. 1]) 

(5.17) \F(x,t)\ < Clxl-T^-M 2 , for all ()<??< 4 

and some constant C > independent on x, t and on < a < 4. Moreover, 

(5.18) F(x, t) = r 2 F(x/Vi, 1). 
These bounds imply the useful estimates 

(5.19) \\F(t)\\ p < Ct~ 2+ ^ (l<p<oo). 

Applying the first of ()5.6[) with p = q = 6, we get ||u(t)||e < i -1 / 4 . Similarily, ||i>(s)||oo < 
(1 + i)~ 1/2 . Hence, 

||B(« J «)|| 0O < C [ \\F(t-s)\\ 6/5 \\u®v(s)\\ e ds 



(5.20) 



< C (\t- s )-^ s -V\i + s )-^ds 
Jo 

< cii + ty 1 / 2 . 



It remains to establish a pointwise estimate in the in the region {(x, t) : \x\ > 2yi}. Let 
us decompose 

B( u ,v) = r 1 + r 2 , 

by splitting the integrals as JoJL^Mrt • • • an d Jq J\ y \>\ x \/2 • • • • For the estimate of I[ we 
use |u| < s -1 / 2 , \v\ < |y| _a s( a_1 )/ 2 and \F(x,t)\ < C\x\~ z t~ l l 2 . For the estimate of I' 2 we 
use again \u\ < s -1 / 2 , \v\ < C|y|~ a ,s( a ~ 1 )/ 2 , and the L 1 -estimate for F (see (|5.19p ). This 
leads to 

\B(u,v)\(x,t) < Ct\x\~ a &- l V 2 . 

The conclusion follows combining this with estimate (|5.20p . The proof of estimate ()5.14p 
is similar. Notice the limitation (2a)* < 4, which is due to the restriction on tj in in- 
equality (]5.17p . The proof of (15. 15f> also follows along the same lines and is left to the 
reader. 

□ 
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We finish with B(8,u). 
Lemma 5.7. Let a > 1, b > 3. For some constant C > depending only on a,b, we have 

(5.21) I|£(MIU <cHUMU 

and 

(5.22) ||i?(0,it)L < clhIL IML • 

Moreover, 

(5.23) \\B(6,u)\\ y < Cllull^Helly . 

Proof. As before, we give details only for the first estimate. Denoting F(x, t) the kernel 
of e* A V, we can write 

(5.24) B(9,u)(x,t) = J J F(x-y,t-s)(0u)(y,s)dyds, 
Notice that F rescales exactly as F. Moreover, 

(5.25) \F(x,t)\ < C v \x\~n~^~^/ 2 , for all < rj < oo 

These are the same estimates as for F, but there is now no limitation to the spatial decay 
rate (i.e., the restriction r\ < 4 appearing in (15. 17ft can be removed). 

Therefore, the space-time pointwise decay estimates for B(9,u)(x,t) can be proved 
essentially in the same way as in the previous Lemma. 

The L^estimate (useful for estimating the y-norm is straightforward: 

< C [\t - a)- 1/2 |K*)||oo||0(a)||i ds < C\\u\\ x \\0\\y. 
Jo 

This allows us to conclude. 

□ 

Proof of Proposition 12.41 We need two elementary estimates on the linear heat equa- 
tion. Namely, 

(5.26) l|e* A 6>ob < c(ll^o||i + ess sup |x| 3 |6» (x) 

* X 

and 

(5.27) ||e* A «o||^ < Cesssup|x| |uo(x)|. 

X 

Both estimates immediately follow from direct computations on the heat kernel gt(x) = 
(47rt)" 3/2 e" |:r ' |2/(4 * ) . (See, e.g. [UES]) Here one only needs to use \g t (x)\ < C\x\~ 3 and the 
usual L 1 -L 00 estimates for g t . Letting U = e tA UQ and = e tA #o> by assumption (|2.10j) we 
get, for some C > 0, \\U\\x + \\@\\y < Ce. The system (|5.ip can be written in the abstract 
form (]5.2p . By inequalities (I5.8P , (I5.15P and (15.23H , all the assumptions of Lemma 15.11 are 
satisfied provided e > is small enough. The conclusion of Proposition 12.41 readily follows. 

□ 
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Proof of Part (a) of Proposition 12.51 By construction, the solution (u, 0) of Proposi- 
tion [23] is obtained as the limit in X x y of the sequence (u k ,9 k ) defined in f|5. 3j) . 

By the first of assumptions (|2.12p . and applying straightforward estimates on the heat 
kernel (see also @1|29]), \e tA u (x)\ < C(l + |x|)" a and |e' A n | < C(l + t)~ a / 2 . (Here we 
need < a < 3). These two conditions imply in particular that e tA uo 6 X a . Similarily 
one deduces from the second inequality in (|2.12p that e* A #o G 34 (when 6 = 3 one here 
needs also 9q £ L 1 ). 

By estimate (|5.2ip and assumption (|2.10p . for all k = 1, 2 . . . we get 

\\9 k+1 \\y b <\\e tA 9 \\y b +Ce\\e k \\y b . 

If e > is small enough then Ce < 1 (the size of the admissible e thus depend on b and, 
as we will see later, also on a). Iterating this inequality shows that the sequence (9 k ) is 
bounded in 3V 

Combining estimates (|5.7p with (|5.13p we get 

\\u k+1 \\x a < \\e tA u \\ Xa +Ce\\u k \\ Xa +C\\9 k \\y,. 

Assuming Ce < 1, we deduce from the boundness of (0 k ) in 34 that (u k ) is bounded in X a . 
Thus, the solution (u, 6) belongs to X a x 34 and the first part of Proposition 12.51 follows. 

Proof of Part (b) of Proposition 12. 5L The proof of the second part of Proposition (|2.5p 
is quite similar but relies on the use of slightly different function spaces. So, let a > 2. We 
define X a as the Banach space of divergence vector fields u = u(x,t) such that, for some 
C > 0, 

(5.28) \u(x,t)\<C inf \x\^{l + t) { ^ 2 ^ 2 . 

0<ri<a 

For b > 4 we define the space 3^ of functons 6 £ L^{L\) satisfying the estimates 

(5.29) \e(x,t)\<C inf |x|-"(l + t)^l 2 . 

0<r)<b 

Such spaces are equipped with their natural norms. 

Notice that the spaces X a and 34 differ from the their counterparts X a and 34 only by 
the fact that the time decay conditions are slightly more stringent in the former case. 

The counterpart of estimates (15. 6|) are 

\\u(t)\\ LP , q < C\\u\\x a {l +ip ( §~ 2) , max{l, §} < p < oo, 

\\o(t)hp'« < c\\e\\y 4 (i + t)^l- 4 \ i < P < oo. 

The first of (|5.30p can be completed by 

(5.31) IK*)||i < C7||u||^ o (l + if3<a<4. 

This last estimate follows immediately by splitting the integral f \u\ into the regions 
\x\ > t x l a and |x| < t 1 ^. 

We also notice the continuous embedding 

(5.32) y b cL?(L\), if 6 > 4 



(5.30) 



\L 
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that follows easily by splitting the integral J* (1 + \x\)\9(x, t)\ dx into L x \<(i +t \i/2 ■ ■ ■ an d 
iix|>(i+t)i/2 an d using the bound \9(x,t)\ < C(l + t)~ 2 for the first term and |#(x,t)| < 
C|x| -6 for the second one. 

Lemma 5.8. Let 2 < a < 4. If f 9(t) dx = for all t, then for some C > 0, 

(5.33) \\L(0)y a <C\\9\\^. 

Assume, without restriction, = 1- The time decay estimate 

/ K(t - s) * 9(s) ds ^Cil + t)- 1 
Jo 00 

immediately follows by the L?-I? Young inequality. It only remains to prove that L{9) 
can be bounded by C|x| -a (l + i)( a ~ 2 )/ 2 for all (x,t) belonging to the parabolic region 
\x\ > 2\ft. Thus, we decompose 

(5.34) L{9) = I 1+ I 2 + I 3 = + I 1)2 + J lj3 ) + J 2 + I 3 . 

where I\ , I 2 and I 3 are as in Lemma 15.51 and and the terms contributing to I\ are defined 
below. First, 

(l lt i).(x,t) = K ji3 (x,t-s) J 9{y,s)dyds = 
by the zero- mean assumption on 9. Next 

(h,2)Ax,t) = - K jj3 (x,t-s) 9(y,s)dyds 

J JO J\v\>\*\/2 

and 

(h,3)j(x,t) = - / (/ VK J>3 (x- Xy,t-s)dy ) ■ y 9(y, s) dXds, 

J0 J0 \J\y\<\x\/2 / 

where we have used the Taylor formula to write the difference ¥Lj -3 (x—y, t—s)—M>j t3 (x, t—s). 
From the bounds \K(x,t)\ < C|rr| -3 and \S/K(x,t)\ < C|x| -4 we get 

\h,2\ + |/l, 3 | < C\ X \-\ 

where we used the continuous embedding of L\ into 3^4 that follows from the definition of 
such space. The above pointwise estimate is even better, in our parabolic region \x 
than what we actually need. 
Next, 

(h)j = / K jt3 (x - y,t - s)9(y,s)dyds 

JO J\x-y\<\x\/2 

can be bounded as follows 

r-t 

-4 



|/2|<C|xr 4 / / \K l3 (x-y,t-s)\dyds 

Jo J\x-y\<\x\/2 

< C\x\~Hlog(\x\/V~t) 

< C\x\- a t( a ~ 2 V 2 

for all (x,t) such that \x\ > 2\/J (recall that 2 < a < 4). 
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Moreover, using again \K(x — y,t — s)\ < C\x — y\ 3 , shows that 

(h)j= [ [ K jj3 (x-y,t-s)9(y,s) 

JoJ\y\>\x\/2, \x-y\>\x\/2 

can be bounded by 

c\x\~H. 

So far, we proved that 

\L(9)\{x,t) < C\x\- a t^ a - 2 ^ 2 < C\x\- a {l + t)^ a - 2 ^ 2 

for all (x,t) such that |x| > 2\ft. Our previous L°°-bound on L{9) implies the validity of 
such estimate in the region \x\ < 2\ft. We thus conclude that L(9) £ X a and Lemma HT81 
follows. 

□ 

Next Lemma is a simple variant of Lemma 15.61 
Lemma 5.9. Let 2 < a 
(5.35) 

and 

(5.36) \\B(u,6)\\y b < C\\u\\ x \\9\y 

Proof. This Lemma can be easily proved following the steps of estimates (|5.13p and (|5.2ip . 
We thus skip the details. 

□ 

The last estimates that we need, concern the heat equation. The computations are 
straightforward (see [U [29] ) . Recall that 2 < a < 4 and a / 3. Moreover, we assumed 
1^0(^)1 < C(l + |x|) _a . When 2 < a < 3 then, as we already observed, \e t/ ^uo(x)\ < 
C(l + |x|)~ a and |e' A uo| ^ C(l + t)~ a / 2 . In fact this estimate remains valid also for 
3 < a < 4 (here one uses that uo is integrable and divergence free, and so J uq = 0). Thus, 
in particular, e tA UQ G X a . 

In the same way, one proves that by our assumptions e* A #o £ f° r b > 4. 

Therefore, going back to the approximation scheme (15.3(1 and arguing as in the proof 
of Part (a) of Proposition (12. 5j) we see that the sequence (9 k ) is bounded in 3^ and (u k ) 
is bounded in X a . Part (b) of Proposition (|2.5p follows. 

□ 



< 4 and b > 4. Then, for some constant C > 0, 
\\B(u,v)\\x a < ClIuH^llull^ 



6. Asymptotic profiles and decay of strong solutions 
We denote by E(x) = -p, the fundamental solution of the Laplacian in M 3 and by 

\x\ 

{E Xj:Xk )(x) its second order derivatives for 1 / 0. Notice that E XxjX3 is a homogeneous 
function of degree —3. Next lemma describes the asymptotic profile for L(9)(x,t) as 
\x\ — > 00, by establishing that 

L(9)(x,t) ~ ( / 6»o ) t(E Xj:X3 )(x), as |x| » Vt. 



DECAY AND GROWTH OF SOLUTIONS OF BOUSSINESQ (July 29, 2010) 27 

Lemma 6.1. Let 9 = 9(x,t) be any function satisfying the pointwise estimates (|2.14[) . for 
some 3 < b < 4 and such that J 9(t) = f 0q for all t > 0. Then the j -component L(9)j(t) 
of L{9) can be decomposed as 

(6.1) L(6)j(x,t) = (J }t(E Xj , X3 )(x) + K'fat) (j = 1,2,3), 
where the remainder function TV satisfies, 

(6.2) \n\x,t)\ < C\x\~ b fl-VI 2 log(|x|M), V (x, t) s.t. \x\ > 
In particular, in the region \x\ 3> \ft, one has 

\K'(x,t)\ ^ Ct\E XjiXk {x)\ 

along almost all directions. 

Remark 6.2. This idea of obtaining informations on the large time behavior of solutions 
by first studying their behavior in the parabolic region \x\ ;§> \/t comes from [6]. 

Proof. We go back to the decomposition (|5.34p of L(9), as done in Lemma 15.81 We now 
treat I2 using the estimate \9(x,t)\ < C\x\~ b (1 + i)^" 3 )/ 2 . This yields to the inequality, 
valid for \x\ > 2\/t, 

\h\(x,t) <C\x\-\l + t)W 2 iog{\x\/y/i). 
With the bound on 9 we obtain also 

IhKx^KClx^il + tf-^ 2 . 
Now recall that l\ = I\ t i + I12 + *i,3> where 

(6.3) (h,3)j(x,t) = - [ [ ([ VK jt3 (x-Xy,t-s)d\) ■ y9{y, s) dy ds. 

Jo JO \J\y\<\x\/2 J 

and that VK satisfies the estimate |VK(x,t)| < C|x|~ 4 . Then, since 3 < b < 4, 
1*1,31 (M) < C\x\- 4 f I \y\ \8(y,s)\dyds 

JoJ\y\<\x\/2 

The estimate for (lip)] = — Jq K?,3(#, t — s) f\ y \>\ x \ ji 0{y-) s ) dy ds, is straightforward: 

|*l, 2 |0M) < C\x\- b {l + ti b -^' 2 . 
Finally, since the mean of 9 remains constant in time, 

(h,l) j (x,t) = J K ji3 (x,t-s) J 9{y,s)dyds = (J 9 ^j J K j>3 (x, t - s) ds. 

But the following decomposition of the kernel K, established in [7], holds : 

K j:k (x,t) = E Xj!Xk (x) + \x\- 3 ^ j:k (x/Vi), j,k = 1,2,3 

where is fast decaying: |^(y)| < Ce~ c ^ 2 for all y G R 3 and some constats c, C > 0. 
Hence, we can estimate |^(y)| < C|y|~ b+3 . 
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We now define lZ'(x,t) through the relation 
L(9)(x,t) ■- 



J 0^\t{E XjiXz ){x)+K'{x,t) 



and all the previous estimates imply \lZ'(x,t)\ < C\x\ b (l + t)( b x ^/ 2 log(|x|/\/^) for all 
\x\ > 2y/t. 

□ 

In the case J 9q = 0, we can use the following variant of Lemma 16. 11 

Lemma 6.3. Let 9 = 9(x,t) be any function satisfying the second of (|2.16p . for some 4 < 
b < 5 and such that J 8(t) = J 8q = for all t > 0. Then the j -component L(9)j(t) of 
L(9) can be decomposed as 

(6.4) L{9)j{x,t) = -VE XjXz {x) ■ Qf y y9(y,s)dyds} + K"j(x,t) 
with 

(6.5) \n"\(x,t) < C\x\- b (l + t) {b - 2)/2 log(\x\/Vi) 
for all (x,t) such that \x\ > 2\ft. 

Proof. We only have to reproduce the proof of the previous Lemma with slight modifica- 
tion. Using the estimate \9(x, t)\ < C|x| _b (1 + t)( fe ~ 4 )/ 2 we now obtain, for \x\ > 2y/t, 

\h\(x, t) < C\x\- b (l + t)(*- 2 )/ 2 log(|x|M), 

next 

\h\(x,t)<C\x\- b (l + t)^ b - 2 ^ 2 . 

and 

\h,2\{x,t) <C\x\- b {l+4 b - 2 )/ 2 . 

Next, by the vanishing mean condition I\ \ = 0. It remains to treat J13. We can 
decompose I\ j3 , whose j-component we recall is 

n[K jt3 (x -y,t- s) - K j}3 (x, t - s)]9(y, s) dy ds, 
J<M/2 

into the sum of three more terms 

^1,3 = ^1,3,1 + -^1,3,2 + h,3,3- 

Such decomposition is performed exactly in the way we did in the proof of Lemma 16.11 
The j-component of the first term is thus 

(6.6) (Ii,3,i)j(x,t) = -J VK jj3 (x,t - s) ■ J y9(y,s)dyds. 
The second term, 

(h,3,2)j(x,t) = / VK jj3 (x,t- s) ■ / y9(y,s)dyds, 
Jo J\y\>\x\/2 
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can be bounded by the right-hand side of (|6,5p using |VK(x, t)\ < C\x\~ 4 and \9(x,t)\ < 
C\x\~ b (l + i)( fe ~ 4 )/ 2 . Next, the j-component of Ji )3j3 , 

n[K jt3 (x -y,t- s) - K j>3 (x, t- s) + VK(x, t - s) ■ y]9(y, s) dy ds, 
,j\<\*\/* 

can be treated with the Taylor formula. The simple estimate |V^IK(x,i)| < C|x| -5 allows 
us to see that also h j3:3 is bounded by the right-hand side of (|6.5|) , Therefore, both -^1,3,2 
and Ii,3,3 can be included into the remainder term lZ"(x,t). 

Let us go back to (|6.6j) . As shown in [7], the following decomposition holds true: 

(6.7) VK j>k (x,t) = VE XjXk {x) + \x\-^(x/Vt), 

with < Ce~ c \ y \ 2 for some constants C, c > and all y G M 3 . In particular we can 

estimate |^(y)| < C|y| - ^ -4 ^. On the other hand, J \y\ \6(y,s)\dy is uniformly bounded 
because of the embedding (|5.32p . This shows that (11,3,1)3(2;, i) can be written in the 
region |x| as in the right-hand side of (|6.5p (even without logarithmig factors). 

This finally gives (16. 4p . 

□ 



We can now establish our main results as simple corollaries: 

Proof of Theorem I2.6L part (a). Let (u, 9) be a mild solution of the system (15. ip . 
satisfying the pointwise decay estimates (I2.13p - (l2.14p . with a > | and b > 3. Recall that 
the spaces X a and 3^ decrease as a and b grow. Without restriction we can then assume 
I < a < 3 and 3 < b < 4 in our calculations. According to our notations, we can write 

(6.8) u(x, t) = e tA u (x) + B(u, u)(x, t) + L(6)(x, t). 

By estimate (|5.14p . owing to the condition a > |, we have 

lim *(";")(*.*) =0 , 

Therefore B(u,u) can be included inside the remainder term in the asymptotic profile of u 
for J^i — > 00. Moreover Lemma (|6.ip and the condition b > 3 guarantee that 

L(9)(x,t) = (J o) tVE X3 (x) + o(t\x\~ 3 ), as ^ -> oo. 

This yields the asymptotic profile (|2.17p for u. 

Let us prove here also the claim in Part (a) of Remark l2.71 As usual, we denote gt(x) = 
(4vrt)- 3 / 2 e-l a: l 2 /( 4 *) the standard gaussian. Under the additional assumption |uq| < C\x\ 3 , 
we have J\ y \>\ x \/29ti x ~ y)\ u o{y)\ dy < C|x| -3 . Moreover, by the third of (|2.10p . we have 
also, e.g., K(x)| < C\x\~ 2 . So f\ y \<\ x \/ 2 9t(x - y)\u (y)\ dy < C\x\ sup\ z _ x] <\ x \ /2 g t (z) < 
\ft |^| -3 - Combining these estimates we get 

\e tA u (x)\ < t|x| -3 , as t > 1. 

Thus, when f 9q ^ 0, under the additional condition, |uo(x)| < C|x| -3 , for |x| » yi >■ 1, 
the solution u behaves like (f uo)tVE X3 (x) along almost all directions. More precisely, 
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using spherical coordinates and letting x = pco, with p > 0, for almost all u in the unit 
sphere, we have 



J^«iOM)/(C/ u o)tE Xj , Xk {x) j = l. 

Proof of Theorem 12. 6L part (b). We can assume without restictrions for the calcula- 
tions below that that 2 < a < 3, and b > 4. Moreover, by our assumption J 8q = 0. The 
velocity field n belongs to A^. Then 

|B(u,«)|(x,t) < C\x\~H 1/2 

as it can be proved easily by splitting B(u,u) into Joi|?/|<|xi/2 an< ^ JoJh/|>|a;|/2 anc ^ usm S 
|n| < C|x| -2 and ||u(i)|| < Ci -1 / 4 (see (|5.30p ). We get easily a bound for both terms in 
the region \x\ > 2y/i, implying 

lim l*(";")IOM) =0 
(M^ m t\x\ 

Thus, B(u, u) can be included inside the remainder term. Applying now Lemma [6.31 yields 
the asymptotic expansion (|2,19p . 

Under the additional assumption | no (x) | < C|x|~ 4 . we get, by (|2.10p . |no(x)| < C|x| -5 / 2 . 
If we use the bound |<7t(x)| < Ct 3 ^\x\^ 9 ^ 2 and the usual L 1 -estimate for gt we get 
|e tA a(x)| < C(l + i 3 / 4 )|x|" 4 . Thus, 

|e* A n (x)| < t|x|~ 4 , as t > 1 

and the last claim (made rigorous exactly as above) of the theorem follows. 

□ 

We now deduce from Theorem 12 . 61 sharp upper and lower bound estimates in L p -spaces. 

Proof of Corollary 12.81 part (a). The upper bounds are simple: indeed, appying the 
arguments that we used in the proof of Lemma 15.41 to (1 + | • |) r |n(-, t)\ instead of n (and 
putting p = q, in a such way that Lorentz spaces boil down to the usual Lebesgue spaces) 
gives the result. 

We now discuss lower bounds. By the proof of Lemma 16. II and of Theorem 12. 6} we can 
find an exponent ij > (any < rj < min{2a — 3, b — 3, 1} will do) such that, for j = 1,2,3, 



\uj(x,t) - e tA (n j)(x)| 



> t 



je \E Xj>X3 (x)\-Ct\x\- 3 {^) 



provided |x| > Ay/i, and A > is large enough. 
Consider the parabolic region 

T> A ,t = {(x,t): \x\ > Ay/t}. For 1 < p < oo, we denote 

by II • h?.(D Ait ) the norm 



mv A , t ) = (j v \f(x)\ p {i + \x\rd x y p 
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Then for all t > 1 and 1 < p < oo, r > such that r + | < 3, 



itj(£) - e iA n ,j 



> Ct 



if(»A,t) 



>Ct^ r+ I^A r+ t 3 (C' 



L P r(V A , t ) 



Ct 



This shows that it is possible to define a continuous function (j>: [0, e 2 ] — >■ M + , (where e > 
is the constant of Proposition 12. 4p such that (f)(0) = 0, <j> is strictly positive outside the 
origin, and satisfying 



(6.9) 



u(t) 



e tA u 



> 



u(t) 



e tA u Q 



><P{\f0o\)t 



1(^+1-1) 



for all t > large enough. By comparing the two terms inside the parentheses in the 
inequality above, i.e., by taking A such that A -7? < ^-| / #o|> we get an explicit behavior 

for <fi(cr) near zero, namely, <J)((t) ~ c£r 1+: ^ 3 r p\ as <r — )• + . with c > small enough. 

We now restrict us to the smaller range < r + - < min(3, a), always with r > and 
1 < p < oo. Let us compute the L^-norm of e' A uo- From |uo(x)| < Cnrrn{|a;| — 1 , |x| _a } 
we obtain, for t > 1, e' A//2 |uo| (x)(l + |x|) r < C(l + |x|)~( a ~ r ). Applying the semigroup 
property of the heat kernel, we get, for t > 2, e* A |«o|(x)(l + |x|) r < e* A//2 (l + |x|)~^ a_r \ 
Computing the L p -norm of this quantity, we deduce 



(6.10) 



|e* A u || L p < \\g t / 2 * (1 + I • I)" 

< C||fl , t/2|U°'.P 
= Ct -±{a-r-l)_ 



-(a-r)i 



with 1 + i 

p 



+ 



In this computation, denotes as usual the weighted LP space, whereas L a ' p is a Lorentz 
space. Here we made use of Young convolution inequality, generalised to Lorentz spaces 
(see [271 Prop. 2.4]). 

By comparing the large time behavior of the RHS in expressions (|6,9|) - (|6.1(J|) . we deduce 
the lower bound 



Ht)\\L*>M\feo\)t-^ r+ * 



for all t > to 



where to > is some constant depending on all the parameters and the initial data, but 
independent on t. 

□ 



Proof of Corollary 12. 8|, Part (b). The estimate from above follows applying inequali- 
ties ff5730^ and fl&E]) to (1 + | • \) r u. 

Let us now estimate ||w(£)||x,p from below. The proof is based on the asymptotic expan- 
sion (|2.19p . Computing the third order derivatives outside the origin of the fundamental 
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solution E of —A in R 3 , i.e., E{x) = Q, shows that (see also Eq. (9b)]) 

r(Y) „ ,..m.. ,2 



( 6 - n ) E x x hiXk (x 



Oj,h,k{x)\x\ z - 5xjX h x k 



j,x h ,x k \-"} ^3/2 \ x \7 



with (Jj,h,k{x) = 5j,hXk + $h,kXj + 6k,j%h- It is now easy to see that the expression 
Y2h k=l Exj,x h ,x k {x)Mh,k identically vanishes if and only if Mh t k is a scalar multiple of 
the identity matrix. Let m(t) = f^J y6(y, s) dy ds. We deduce that the homogeneous 
function of degree —4 

V/;, ,.. r; (.r) •.,,( /) 
is identically zero, for any fixed t, if and only if |m(i)| = 0. 



l ft 
t 



Jol yG{y,s)dy ds 



7^ 0, then there exists c > such that, for all t 



•Mr+i-2) 



If m = lim inf t-^oo 
sufficiently large, 

|| V^,,^ • m(t)\\ L P (VAt) > cmt^r 

The condition on the remainder 7Z obtained in Theorem 12. 61 then implies, some constant 
c > and for all t large enough, 

ll^-e^olL^^^+f- 2 ) 

with r>0, l<p<oo tels que r+- < 4. It remains to prove that, when r+- < min{a,4} 
then we have 

(6.12) l|e* A n || L p = o(^ (r+ t~ 2) ), for t -)■ oo. 

Recall that we assumed a > 2. When 2 < a < 3, we can simply use inequality (|6.10p . 
When a = 3 there is nothing to prove because we reduce to the previous case by picking a' , 
with r + 1 < a' < 3. So, consider now 3 < a < 4. In this case uo is integrable and f uq = 
by the divergence- free condition. Then, for t > 1, e* A/2 |u |(2:)(l + |x|) r < C(l + |x|)-( a - r ). 
Thus, 



5-13) \\e u \\ L p < \\g t /2 * (1 + 



When a — r < 3 we can apply Young inequality in the same way as before and still obtain 
estimate (|6.10p . When a — r > 3, the above quantity is bounded by Ct 2 k p ' (and by 
C v t~ 2{ p~ 3+v) for all r) > when a-3 = 3). In any case, (|6,12p holds true. This establishes 
estimates (j2.24[) 

□ 

Theorem 12.61 has another interesting consequence, that clarifies the importance of the 
restriction r + 3/p < 3 in our previous statements. 

Corollary 6.4. Let (6,u) be a solution as in Part (a) of Theorem \2.6\ We assume, in 
addition, that J 9q 7^ and that the initial velocity satisfies \uq{x)\ < C(l + |x|) _a , for 
some a > 3. Then for all 

3 

r > 0, 1 < p < 00, r H — > 3 

P 
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and for all t > we have 

(6.14) \\u(t)\\ L? =oo. 

Proof. Indeed, E Xj)X3 (x) is a homogeneous function of degree 3, smooth outside the origin. 
Then we can find an open conic set r C R 3 such that \E XjtX3 (x)\ > C|:c|~ 3 > for all 
i £ T, i / 0, for some C > 0. Indeed, by Eq. (l2~T7j) . for \x\ > Ay/i with A > large 
enough and x £ T, and t > 0, we have 

C 

\uj(x,t)\ > —t 
This implies that Uj has an infinite || • ||x£(D A nr) -norm. 

□ 




7. Additional remarks and comments 

In this section we collect a few technical remarks on the main results. These are essen- 
tially small variants of our statements that can be easily proved with minor modifications 
to the proofs. 

Remark 7.1 (on Theorem 12. 2\i . Part (b) of Theorem 12.21 can be streghened as follows. 
Under the same assumptions on 9$ an d replacing the assumption no £ E\ n L 3//2 with the 
weaker condition ||e* A no|| 2 < C(l + t)~ s , for some s > 0, (when no £ L^nL 3 / 2 , this holds 
with s = 1/2), we have 

(7.1) ||n(t)|| < C(l + t)- s *, with s* = min(s, 1/2). 

The above decay condition on e* A no could also be restated in terms of Besov spaces. 

Remark 7.2 (on Theorem 12. 6p . The second term in the RHS of (|2.19p is bounded by 
C|x| _4 t. If one is interested in studying the asymptotic behavior of n only as j^j — > oo 
(with t > not necessarily large), then an additional term in the RHS of (|2.19p should be 
added: 

(7.2) -\7E Xh)Xk : J J (u h u k )(y,s)dyds 

(the : notation means that the Ylhk s y mD °l has been omitted). Such term is bounded 
by C|x| _4 t 1//2 , thus justifying its inclusion inside the remainder when \x\ ^> \/t ^> 
1. With this additional term, the condition of the remainder can be simplified into 

v t 

Remark 7.3 (on Corollary 12. 8p . The restriction r + | < a in Part (a) of Corollary 12.81 is 
natural beacause the decay assumption on no guarantees that e* A n G L v r exactly for those 
r, p satisfying such restriction. However, estimates (|2.23p remain valid in the whole range 
< r + | < 3 (and under the conditions of Part (b) even for 0<r<|<4ifwe want to 
estimate \\u(t) — e t ^uo\\ L p instead of ||n(i)||£P in that expression. 
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